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SYNTHETIC-APERTURE  SONAR:  PERFORMANCE  ANALYSIS  OF 
BEAMFORMING  AND  SYSTEM  DESIGN 


by 

Ezio  G.  Pusone  and  Lewis  J.  Lloyd 


ABSTRACT 


Studies  undertaken  to  assess  the  feasibility  of  applying  synthetic-aperture 
techniques  to  sonar  have  been  based  on  the  theories  developed  for  radar 
systems.  They  have  been  limited  primarily  to  short-range,  high-frequency 
active  systems  because  the  ratio  of  propagation  speed  to  platform  speed  is 
much  smaller  than  that  encountered  with  radar.  This  study  places  emphasis 
on  the  beamforming  aspects,  the  most  important  part  of  a  synthetic-aperture 
system.  A  general  approach  is  followed  to  develop  a  computer  model  to  give 
estimates  of  beam  patterns  with  active  and  passive  modes  of  operation  for  a 
stationary  target.  The  model  allows  the  investigations  of  a  number  of 
features  of  the  performance,  for  example  the  effect  of  random  phase  and 
amplitude  fluctuations,  of  signal  coherence  time  in  the  passive  mode,  and 
of  both  deterministic  and  random  positional  errors  of  the  synthesizing 
array.  In  addition,  the  usefulness  of  weighting  factors  to  steer  and  shape 
the  synthesized  beam  is  considered.  The  results  presented  give  theoretical 
guidelines  for  the  design  of  active  and  passive  synthetic-aperture  sonar 
systems. 
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INTRODUCTION 


The  aim  of  the  synthetic-aperture  technique  is  to  obtain  a  high  angular 
resolution  from  a  small  moving  physical  array  by  combining  coherently  the 
data  from  successive  positions  of  the  small  array.  It  has  been  used 
successfully  in  aircraft  and  satellite  active  radar  systems  for  which  the 
speed  of  signal  propagation  is  very  much  higher  than  the  platform  speed  (1, 
2). 

In  applications  to  sonar,  the  ratio  of  propagation  and  platform  speeds  is 
very  much  lower;  consequently  there  are  sampling  problems.  To  date  the 
applications  in  sonar  have  concentrated  on  uses  for  sea-bottom  mapping  (3, 
4,  5)  and  mine-hunting  systems  [6], 

Previous  studies  on  applications  to  sonar  [7,  8,  9,  10,  11,  12]  have  tended 
to  take  directly  the  results  developed  for  radar,  and  have  concentrated  on 
short-range,  high-frequency  active  systems.  The  present  study  builds  up, 
for  longer  ranges  and  low-frequency  systems,  the  beamforming  process  from  a 
straightforward  consideration  of  the  phase  of  the  received  signals  at 
successive  positions  of  the  physical  array.  This  very  basic  approach  is 
used  to  develop  a  computer  model  on  SACLANTCEN' s  UN1VAC  1100/60  that  allows 
various  aspects  of  synthetic-aperture  beamforming  process  to  be 
investigated  for  a  narrow-band  passive  system  and  for  a  multi-ping  active 
sonar.  The  target/source  is  assumed  to  be  stationary  in  both  cases. 

The  computer  model  allows  the  investigations  of  the  effects  of  the 
following  features: 

a.  Beam  pattern  of  the  physical  array. 

b.  Complex  weighting  functions  for  shaping  and  steering  the 
synthesized  beam. 

c.  Lack  of  coherence  in  the  received  signals  (passive  mode  only). 

d.  Random  perturbations  in  the  amplitude  and  phase  of  the  received 
signals. 

e.  Placement  errors  due  to  a  non-linear  course  of  the  physical  array 
and  to  random  positional  errors. 

The  lack  of  coherence  in  the  passive  mode  of  operation  referred  to  in  (c) 
relates  to  a  temporal  (and  possibly  spatial)  loss  of  coherence  as  the 
synthesizing  array  progresses  along  its  synthesizing  sequence.  This  is  an 
area  for  which  no  published  data  in  the  open  literature  are  available. 
Consequently  a  fairly  simple  assumption  has  been  used  in  this  report, 
namely  that  the  coherence  falls  off  exponentially  with  time. 

The  amplitude  and  phase  fluctuations  referred  to  in  (d)  could  be  caused  by 
a  number  of  factors,  such  as  changes  in  the  reflecting  characteristics  of 
the  target  as  the  synthesizing  sequence  progresses,  fluctuations  in  the 
propagation  path,  or  the  presence  of  noise.  In  order  to  quantify  these 
effects  they  have  been  related  to  the  signal-to-noise  power  ratio  at  the 
output  of  the  synthesizing  array. 
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Two  types  of  placement  errors  (errors  between  the  actual  and  assumed 
positions  of  the  moving  synthesizing  array)  have  been  investigated.  One  is 
a  purely  random  positional  error  caused  by,  say,  random  navigation  errors, 
and  the  other  is  a  deterministic  error  caused  by  the  array  following  a 
non-linear  course  that  oscillates  about  the  assumed  straight-line  track. 
This  is  the  type  of  situation  that  can  easily  occur  at  sea  when  a  ship 
attempts  to  follow  a  straight-line  course  but  invariably  "oscillates" 
around  it. 

We  start  with  an  analysis  of  the  passive  mode  of  operation  because  it  is 
simpler  to  understand;  we  then  consider  the  active  mode. 
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1  PASSIVE  SYNTHETIC-APERTURE  PROCESSING 

In  the  passive  mode  of  operation,  samples  of  the  signal  received  at 
successive  positions  of  the  array  as  it  moves  through  its  synthesizing 
sequence  are  phase-compensated  for  the  delay  time  and  the  movement  of  the 
array  between  samples,  and  then  added  coherently.  The  achieved  performance 
obviously  depends  on  the  phase  stability  of  each  signal  sample  used  in  the 
coherent  addition.  We  consequently  attempt  to  analyze  and  assess  the 
various  possible  sources  of  phase  (and  amplitude)  errors. 

The  main  assumptions  for  the  passive  mode  of  operation  are  that  the 
radiating  source  a)  is  stationary,  b)  exhibits  perfect  spatial  and 
temporal  coherence  that  exists  during  the  synthesizing  period,  and  c) 
carries  out  own-doppler  nullification  before  the  synthetic  beamforming 
process.  With  regard  to  the  complex  weighting  factors  (delay  and  steering 
vectors)  used  for  beamforming,  these  are  computed  on  the  basis  that  the 
synthesizing  array  moves  along  an  idealized  straight-line  course  at  a 
constant  speed  (U)  and  that  the  received  signals  are  sampled  at  a  constant 
sampling  interval  (At). 

On  the  basis  of  these  assumptions,  a  computer  model  has  been  developed 
using  the  equations  derived  in  App.  A.  This  model  allows  the  computation 
of  the  synthesized  beam  pattern  for  the  situations  given  below.  The 
authors  used  the  procedure  described  in  [13], 


1.  1  Idealized  Beamforminq 


Ihe  first  situation  assumed  completely  idealized  parameters,  full  signal 
coherence  for  the  total  synthesizing  period,  no  phase  or  amplitude  errors, 
and  no  placement  errors  (Fig.  1).  The  general  equation  for  the  beampower 
output  is  Eq.  A. 8  of  App.  A: 


r  n 

? 

~N 

p  = 

[k*/k  Bk  cos<»Sk-8Fk) 

+  j 

k2=1\  Bk  sin(0Sk~  W 

In  the  far-field  of  the  synthesized  aperture  this  equation  becomes  the 
standard  beam  pattern  expression  for  a  linear  array  (see  Eq.  A. 14  of 
App.  A). 
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1. 2  Effect  of  Signal  Coherence 

The  second  situation  assumed  simply  and  plausibly  that  the  signal  coherence 
falls  off  exponentially  with  time.  The  expression  for  the  synthesized  beam 
pattern  is  now  given  by  Eq.  A. 18  of  App.  A: 

pVi  s2k  +  2i=lll  w’iexp(‘ltk+i-tkl/Tc)  •  (Ec*' 2) 

In  this  equation  the  parameter  T  is  referred  to  as  the  "coherence  time"  of 
the  signal .  c 

1. 3  Effect  of  Random  Amplitude  and  Phase  Fluctuations 

Random  amplitude  and  phase  fluctuations  could  be  caused  by,  for  example, 
changes  in  the  reflecting  characteristics  of  the  target/scatterer,  by 
propagation  effects,  varying  multipath  structures,  or  the  presence  of 
noise.  Characterizing  these  fluctuations  is  not  easy,  particularly  because 
there  exist  no  published  data  on  fluctuations  for  a  moving  platform. 
Consequently  we  assume  here  that  the  fluctuations  are  a  first-order  Markov 
process.  In  order  to  quantify  the  fluctuations ,  we  further  assume  that 
they  have  a  root  mean  square  (rms)  value  inversely  proportional  to  the 
s igna I- to-noi se  power  ratio  at  the  "hydrophone"  level,  a  mean  value  of 
zero,  and  a  correlation  time  that  varies  exponentially  from  sample  to 


sample  (that  is,  with  time).  Thus,  for  example, 
characterized  by: 

the  phase  fluctuations  are 

~k  =0 

(Eq.  3) 

T2  =(S/N)‘1 

(Eq.  4) 

E{e.*£.A.J  =  £,  •  .  exp(- |t,  .  .-t.  1  /Tc)  , 

*  k  k+r  k  k+i  j  k+i  k|  t  * 

(Eq.  5) 

where  T^  =  correlation  time,  which  is  preselected. 

When  Tg  is  small,  that  is,  when  the  correlation  between  errors  is  small, 

their  effect  is  small.  Alternatively,  when  Tf  is  large  their  effect 
is  greater. 

The  resulting  analytical  equation  for  the  effect  of  random  amplitude  (q) 
and  of  phase  errors  (e)  on  the  beam  pattern  is  given  by  Eq.  A. 34  of  App.  A; 

N  _ 

P  =  2  A*  f2(ak)f2(Rk)(b*+n2)  + 

*v;ii;iAkf(.k)f(Rk).Ak+1f(.  k>j>-ib2  + 

+T2  exp(-|  t  k+j-t  J/T  £  )}  •  { 1-^  exp(  -  )t  k+f  t  |  /\  )}  -cosle^.  )] 


(Eq.  6) 
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1.4  Effect  of  Placement  Errors 

By  placement  errors  we  mean  the  difference  between  the  assumed  and  actual 
position  of  the  synthesizing  sequence;  two  types  are  investigated  (see 
Fig.  2).  first,  the  effect  of  the  placement  error  due  to  random  positional 
errors  of  the  synthesizing  array,  (those  for  example  caused  by  random 
navigational  errors)  is  investigated.  Second,  the  effect  of  the 

deterministic  placement  error  due  to  the  synthesizing  array  following  a 
non-linear  course  that  oscillates  around  an  assumed  straight-line  course  is 
examined.  Because  these  errors  are  in  terms  of  quantities  that  are  not 
dimensionless  their  effect  will  be  lessened  by  reducing  the  frequency  of 
operation. 


FI  .  .  PASSIVE  SAS  -  REASONS  FOR  LOSS  IN  PERFORMANCE 


1.4.1  Random  Placement  Errors 

The  random  placement  error  is  defined  as  a  circular  gaussian  distribution 
of  array  positions  about  the  assumed  position.  The  errors  are  independent 
of  each  other  and  of  the  actual  signal  values.  Thus  the  effect  of  these 
signals  is  independent  of  signal  direction  and  can  be  characterized  by  a 
single  standard  deviation  a(y).  This  can  be  converted  into  an  effective 
variance  of  phase  as  follows: 

VAR  (phase)  =  (u(y)-2nA)2  . 


Consequently,  as  shown  in  Eq.  A36  of  App.  A,  the  effect  on  P  can  be 
obtained  from 

P  =  2  AkBV  2  2  2  AkBkAk+iBk+i(1‘VAR(Phase))-  cos(®k~9k+i }  ■ 

K  1  1  1  K  1  /  •  _  -i  \ 

(tq.  7) 


7 


SACLANTCEN  SR-91 


1.4.2  Effect  of  Deterministic  Placement  Errors 

As  outlined  earlier,  these  are  positional  errors,  and  consequently  phase 
errors,  that  are  caused  by  the  synthesizing  array  following  an  oscillating 
course  around  an  assumed  straight-line  course. 

A  phase  error  is  introduced  (see  App.  A,  Sect.  A4.2)  because  the  actual 
ranges  for  the  steering  and  delay  vectors  are  different  from  those  that 
were  assumed.  I  he  result  is  that  the  ranges  Rj...  and  R^k  defined  in  App.  A 

must  be  modified  by  the  across-track  and  along- track  displacements  of  the 
synthesizing  array  from  the  assumed  positions.  The  resulting  equation  for 
the  synthesizing  beam  pattern,  excluding  all  other  errors,  is  given  by 
Eq.  A. 44  of  App.  A. 


N 

p=  2  a*  b-; 

k=l  K 


}  !AkBkAk+i8k.i  cos(e;-ek+i)  • 

i-l  k=i 


where  0^  ,  ii’  .  are  the  modified  values  obtained  using  tqs. 
and  A.  43  of  App.  ■' 


(Eq.  8) 
A. 41,  A. 42 


1.5  Resull_s 

The  results  in  f  igs.  3  to  7  cover  only  the  main  lobe  and  nearby  sidelobes 
of  the  synthesized  beam.  They  are  in  a  standardized  form,  independent  of 
frequency,  for  a  synthesized  aperture  (L)  of  50X  with  signal  samples  taken 
at  0.4A  spacing  (no  diffraction  lobes  in  real  space  for  all  angles  of 
steer)  and  for  an  array  speed  along-the-track  of  U  =  2  m/s.  The  sine  of 
the  steer  angle  covers  the  interval  ±0.1  and  the  sine  of  the  3  dB 
beamwidth  is  approximately  ±  0.01. 

To  relate  the  curves  to  a  frequency  we  simply  use  the  required  value  of  A. 
to  obtain  L  and  consequently  the  value  of  U  to  obtain  and  the  sampling 

interval  (at  0.4A/U).  Thus  a  speed  other  than  2  m/s  alters  only  the 

values  of  Ts  and  at. 

For  a  synthesized  aperture  other  than  50A  we  alter  appropriately  the 
coordinate  of  the  sine  of  the  steer  angle  to  cover  the  new  range  of  angles, 
and  adjust  the  vertical  coordinate  for  the  different  value  of  gain.  For 
example,  with  a  10\  array  the  coordinate  of  the  sine  of  the  steer  angle 
becomes  ±  0.5  and  the  vertical  axis  is  reduced  by  10  log  5.  Additionally, 
the  value  of  T  is  reduced. 


1.5.1  Effect  of  Coherence  Time 

Figure  3  shows  beam  patterns  plotted  for  different  ratios  of  T  /T  where 
T^  is  the  synthesizing  time.  It  is  seen  that,  as  modelled  in  this  study,  a 
value  of  Tc  at  least  equal  to  or  greater  than  T^  is  necessary  to  retain 
a  reasonable  beam-pattern.  Alternatively,  Ts  must  be  less  than  Tc  . 
Reported  values  of  T£  applicable  to  synthetic-aperture  studies  are  difficult 


8 


SAf.L  ANTCEN  SR  -  9 1 


to  find  in  the  literature;  one  appropriate  set  |14|  toi  a  400  Hz  signal 
spans  values  of  148  s  for  a  range  of  93  n.mi  and  450  s  fur  a  range  of 
268  n.mi.  Some  propagation  measurements  under  the  ice  in  the  Arctic  Jib) 
suggest  that  long  coherence  times  can  exist  there. 


1.5.2  Effect  of  Random  Phase  and  Amplitude  Errors 

Eigure  4  shows  the  results  obtained  by  modelling  random  phase  and  amplitude 
errors  as  a  function  of  the  signal-to-noi se  power  ratio  (S/N)  at  the  input 
to  the  beamforming  process.  Although  the  definitions  of  f,  and  n  are 
such  that  they  can  only  show  trends,  they  are  not  of  an  unreasonable 
magnitude.  With  these  definitions,  the  results  suggest  that  a  S/N  of  at 
least  3  dB  is  required  to  maintain  an  acceptable  muuilobe-to-sidelobe 
ratio.  Alternatively,  phase  fluctuations  of  0.11  of  a  cycle  and  amplitude 
fluctuations  of  0.7  of  the  maximum  (A)  can  be  tolerated. 


1.5.3  Effect  of  Deterministic  Placement  Errors 

Figure  5  shows  the  effect  of  having  a  large  take-off  angle  (y  =  2°)  for 

different  ratios  of  T  /T  where  T  is  the  period  of  use  illation  of  the 
so  o 

array.  A  large  ratio  is  preferable;  a  value  of  5  i  about  the  lowest 
limit.  Obviously  the  value  of  T  depends  very  much  on  i.uw  the  synthesizing 
array  is  moved  through  the  sea  %d  consequently  can  be  defined  only  for  a 
particular  system. 

When  the  ratio  Ts/Tq  is  less  than  or  equal  to  0  i  an  undistorted  beam 

pattern  is  achievable  but  now  a  deflection  of  the  main  lobe  occurs  that 
approaches  y  as  lrs/lQ  is  made  smaller. 

1.5.4  Effect  of  Random  Placement  Errors 

As  seen  from  Fig.  6  and  as  appreciated  from  the  ' pi; vs  its 1  ,  a  much  larger 
random  placement  error  can  be  tolerated  as  the  frequency  is  reduced.  Since, 
in  general,  random  placement  errors  defined  by  a  gaussian  distribution  with 
a  standard  deviation  in  excess  of,  say,  3  m  are  unlikely,  the  effect  of 
this  type  of  error  could  be  made  negligible  by  working  at  a  sufficiently 
low  frequency. 


1.5.5  Effect  of  Weighting 

Two  forms  of  array  weighting  have  been  evaluated  to  show  the  effect  when 
deterministic  errors  are  present.  The  results  are  shown  in  Fig.  7,  in 
which  the  condition  without  weighting  is  for  y  =  2°  and  for  the  worst-case 
ratio  of  Ts/Tq  =  1.  Note  that  the  Hamming  weighting  (controlled  roll-off 

sidelobe  levels)  [161  gives  a  better  result  than  the  Dolph-Chebyshev 
weighting  (constant  value  sidelobe  levels)  (17)  in  that  the  two  peaks  in 
the  beam  pattern  are  eliminated. 
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2  ACTIVE  SYNTHETIC-APERTURE  PROCESSING 

The  system  envisaged  for  active  operation  is  one  in  whu.h  pulse  trans¬ 
missions  are  made  successively  as  the  synthesizing  array  is  moved  forward 
(see  Fig.  8). 

Transmission  is  made  at  periods  corresponding  to  the  assumed  position  of 
the  target  (stationary)  with  own-doppler  nullification  being  carried  out 
before  the  synthetic  beamforming  process.  Consecpient  ly ,  the  periods 
between  the  spatial  samples  are  not  equal,  although  lor  far-field 
beamforming  the  difference  is  very  small.  The  time  between  transmission  is 
dependent  on  the  maximum  range  of  interest  but  is  usually  severely 
constrained  by  the  spatial  sampling  conditions.  Transmission  and  reception 
must  be  made  at  least  at  every  half-wavelength  travelled  in  order  to  avoid 
diffraction  lobes.  This  sampling  distance  can  be  increased  (in  practical 
cases)  to  almost  a  wavelength,  since  in  general  only  small  angles  of  steer 
are  used.  Thus,  for  these  conditions  the  maximum  range  (R)  m  operation  is 
given  by: 


R  =  i  (^) 

max  2  v  U 


(Eg.  9) 


For  a  frequency  of  300  Hz  and  a  platform  speed  U  c.l  2.1  m  =  the  maximum 
range  of  operation  is  approximately  3  km.  This  is  »  c.c  ;  u, libation,  but 
can  be  extended  by  using  a  physical  array  beam  pi!'.  .  nulls  on  the 

direction  of  the  likely  diffraction  lobes  caused  ,  -r  -  imp  liny  when 

using  a  longer  range.  Alternatively,  a  lower  i.s  ,.u  increases  the 
operating  range. 

The  model  also  allows  for  specific  complex  weighting  functions  to  be 
introduced  to  shape  and  steer  the  synthesized  beam.  Ihe  synthesized  beam 
patterns  are  calculated  for  an  ideal  straight  line  course  of  the 
synthesizing  array  and  also  for  the  more  realistic  can:  of  a  non-linear 
course.  A  sinusoidal  course  is  actually  considered  in  the  model.  The 
effect  of  random  positional  errors  on  the  beamfarming  process  is  also 
considered.  In  addition,  an  attempt  is  made  to  include  the  effect  of 
random  fluctuations  in  the  amplitude  and  in  the  phase  of  the  received 
echoes.  The  procedure  followed  is  similar  to  that  used  by  the  present 
authors  in  a  first  study  on  synthetic  aperture  [18J. 


2. 1  Active  Idealized  Beamforming 

The  main  assumption  for  the  active  mode  of  ideal  operation  is  that  the 
phase  of  the  received  echoes  can  be  derived  from  geometrical 
considerations.  This  is  equivalent  to  saying  that  the  received  echoes  are 
not  distorted  on  reflection  from  the  target  by  multipaths,  by  phase 
fluctuations  during  propagation,  or  by  a  low  signal-to-nu i se  ratio. 

The  geometry  of  the  synthesizing  sequence  is  shown  in  Fig.  8.  The 
synthesized  beam  will  be  steered  and  focussed  to  the  target,  which  is 
assumed  stationary.  The  mathematical  derivation  of  active  beamforming  is 
developed  in  App.  B. 


SRBKDUO  MOB  BUMC-NOT  nil MD 

i _  1 
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To  form  a 

weighted  by 
App.  B): 


synthesized  beam  we  sum  N  samples  of  the  received  echoes 
the  appropriate  delay  and  steering  vectors  of  component  dk  (see 


N 

V(P  )  =  2. 
T  k=l 


Vdk‘Wk 


(Eg.  10) 


where  is  the  amplitude  weighting  function  to  shape  the  beam,  or  in  the 

complex  form: 


N 

V(P  )  =  IS 
T  k=l 


(Eg.  11) 


By  varying  the  bearing  fU  of  T  it  will  be  possible  to  trace  out 
the  synthesizing  beam  pattern.  An  estimate  of  the  resulting  beam-power 
output  in  function  of  By  and  f°r  a  fixed  p  ^  is  given  by: 


P(PT>  =  E{V(Pt)-V*(Pt)| 


(Eg.  12) 


Then  the  beam  power  is: 


’  N 

2 

N 

P(By)  = 

_k=l  e  _ 

+ 

*  1(5.) 

.  f  m  K 

Lk  1  J 

(Eg.  13) 


In  App.  B  we  show  that  when  the  ratio  of  Rq  to  the  synthesized  aperture 

length  is  much  larger  than  unity  and  with  uniform  amplitude  weighting  of 
the  sampled  echoes,  the  power  beam  pattern  is  given  by  the  far-field  beam 
pattern  of  the  familiar  form  (sin(Nx)/sin(x))2  [19).  We  also  show  that  the 
3-dB  beamwidth  of  the  active  synthetic-aperture  array  is  one-half  that  of 
the  passive  array  of  the  same  aperture.  This  important  result  can  be 
explained  from  the  mode  of  operation  of  synthetic  aperture.  In  fact,  the 
essential  feature  of  the  synthetic-aperture  technigue  is  that  signals  are 
transmitted  at  successive  positions  of  the  array  and  each  corresponding 
echo  is  delayed  by  a  two-way  path,  and  that  the  corresponding  phase-shift 
is  also  twice  that  experienced  in  normal  array  operation.  Conseguently , 
after  processing  of  the  echoes  the  resulting  beam  pattern  has  a  3  dB 
beamwidth,  which  is  half  of  that  of  a  physical  array  of  the  same  effective 
total  length  (Fig.  9). 
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FIG.  8  SAS  -  ACTIVE  MODE 


9  REDUCTION  OF  UNDERSAMPLING  SIDELOBES  BY  USING 
THE  BEAM  PATTERN  OF  THE  PHYSICAL  ARRAY 
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To  calculate  P(pj)  we  take  the  amplitude  weighting  function: 

W(k)  =  8(k)-b(k)  ,  (Eg.  14) 


where 


B(k)  =  f(ork) -f  (rk) 


(Eg.  15) 


Ihe  function  f(«k)  is  the  amplitude  response  pattern  of  the  synthesizing 
array.  In  our  computer  model  we  have  taken: 


f(u 


k 


where 

and 


sin(nx 

sin(xk 


n 


) 

number  of  hydrophones 

si n  ii.  . 
k 


with 


(Eg.  16) 


(Eg.  17) 


d  =  separation  between  transducers 

A.  =  wavelength 

a  =  angle  (from  broadside)  at  which  the  target  is  seen  by  the 
physical  array. 


The  coefficient  b(k)  is  an  amplitude  weighting  factor  that  can  be  used  to 
shape  the  synthesized  beam  pattern.  Our  model  uses  a  weighting  "Hamming" 
function 


2. 2  Effect  of  Random  Amplitude  &  Phase  Errors 

The  effect  of  random  amplitude  and  phase  errors  on  active  beamforming  are 
evaluated  under  the  same  assumptions  given  in  Sect.  1.3  for  the  case  of 
passive  beamforming.  The  beam  pattern  synthesized  by  active  mode  of 
operation  is  calculated  with  the  amplitude  term  W(k)  and  phase  term  of  S. 
given  by: 

W(k)  =  W(k)  +  g(k)  (Eg.  18) 

?(k)  =  0kT  +  ekT  ,  (Eg.  19) 


where 

W(k)  =  amplitude  with  no  fluctuations 
g(k)  =  random  fluctuation  in  the  amplitude 
0kT  =  phase  with  no  fluctuations 

=  random  fluctuation  in  the  phase. 

For  g(k)  and  ekj  we  have  assumed  the  same  type  of  random  process 
described  in  Ch.  1  for  passive  sonar. 
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2. 3  Effect  of  Placement  Errors 

2.3.1  Effect  of  Deterministic  Placement  Errors 

To  analyze  the  deterministic  placement  errors  in  active  beamforming  we 
consider  the  same  oscillatory  course  of  the  synthesizing  element  as  for  the 
passive  mode  of  operation  (see  Ch.  1). 

The  path  geometry  in  the  active  case  is  shown  in  Fig.  10.  The  error 
produced  in  the  phase  0^  of  the  received  echo  is  given  by: 

2nf 

V  — °  ('■'kT+r"kT)-(RkT+RkT>  '  (Eq'  20) 

kT 

where 

f  =  frequency  of  the  signal 

c  =  speed  of  sound 

r1  ,  r"  =  distance  of  the  oscillating  synthesizing  element  from  the 

target 

R ‘  ,  R"  =  distances  of  the  synthesizing  element  from  the  target  in 

the  case  of  straight-line  course 

Based  on  the  geometry  of  Fig.  10,  we  calculate  the  distances  r1  and  r"  in 
terms  of  the  same  oscillatory-course  parameters  (i.e.  speed  along  the 
course,  maximum  off-track  angle,  period  of  oscillation)  used  in  the  passive 
case.  A  detailed  mathematical  analysis  of  the  phase  error  in  terms 

of  the  path  geometry  described  above  is  given  in  App.  B. 

The  synthesized  beam  pattern  in  presence  of  placement  errors  is  calculated 
from  Eq.  12,  in  which  the  phase  GkT  is  replaced  by  the  sum: 

®kT  =  0kT  +  ^kT  ’  (Eq’  21) 

where 

6kT  =  ideal  phase  term  for  straight-line  course 

ii>kj  =  phase  error  term  given  by  Eq.  20. 

2.3.2  Effect  of  Random  Placement  Errors 

Superimposed  on  the  deterministic  errors  we  have  assumed  that  other 
placement  errors  are  caused  by  random  fluctuations  in  the  motion  of  the 
synthesizing  array.  For  instance,  a  cause  of  random  error  can  be  the 
random  fluctuation  of  the  velocity  of  the  synthesizing  array.  These  random 
placement  errors  have  been  characterized  by  a  circular  gaussian 
distribution  so  that  the  standard  deviation  is  the  same  for  all  signal 
directions.  The  variance  of  the  corresponding  phase  error  of  the  signal  is 
given  by: 


where  \  is  the  wavelength  of  the  signal  and  o  is  the  standard 
deviation  of  errors  distribution.  ^ 
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2. 4  Simulation  Results 

Figure  11  shows  the  result  of  an  active  synthetic-aperture  technique 
working  at  100  Hz  in  the  neat-field.  The  range  R  750  m  implies  a 
sampling  interval  of  about  At  ~  2R/C  ~  1  s.  Taking  U  -  2  m/s  for  the 
speed  of  the  synthesizing  array,  the  corresponding  spacing  d  is  equal  to 
UAt  '2  m.  This  spacing  of  0.4A.  at  300  Hz  is  adequate  to  avoid  grating 
lobes  in  the  total  180°  arc.  The  resulting  synthetic  aperture  has  a  total 
effective  length  of  L  =  T^U  =  400  m.  The  beamwidth  of  the  passive  array  of 

the  same  length  is  given  by 


:-(3  dB) 


b 

~  U  "  '4CT0 


radians  x  57.3 


degrees 
radian "" 


-0.7 ' 


However,  the  3-dB  beamwidth  computed  from  the  cur..,  in  i  ig  11  is  about 
0.35°.  This  result  shows  an  important  feature  of  tin:  ,ynt  lu-tic-aperture 
technique,  namely  that,  the  synthetic-aperture  beaviwidih  is  half  the 
beamwidth  (see  Sect.  2.1)  of  a  passive  linear  array  of  the  same  length. 
This  result  also  differs  from  what  we  would  obtain  f  cm  a  "normal"  active 
array,  which  has  a  beamwidth  that  is  about  forty  pe; .  .  •  (40*)  greater  than 
that  achieved  by  the  active  synthetic-aperture  array. 

Figure  12  shows  that  an  adequate  beam  pattern  can  be  obtained  at  a  longer 
range  of  15  km  at  a  lower  frequency  (f  =  100  Hz)  samples  are  taken  at 
intervals  of  At  ~  20  s  over  a  total  synthesizing  time  of  I  600  s.  An 

effective  synthetic  aperture  of  1200  m,  with  11  -  2  ,n/'s ,  is  therefore 
generated.  The  3  dB  beamwidth  is  0.3°. 

Figures  13  and  14  show  the  effects  of  'deterministic'  placement.  errors  on 
active  beamforming  for  300  Hz  and  100  Hz  operating  frequencies 
respectively.  At  300  Hz  there  is  a  broadening  effect  on  the  main  lobe;  at 
100  Hz  a  double-peak  effect  occurs.  It  appeals  that  the  effect  of 
deterministic  placement  errors  is  more  severe  on  act.,;,  beam;  r.rming  than  on 
passive,  since  the  active  performance  is  acceptable  u.r  values  of  y  <  0.5° 
whereas  values  of  y  <  2°  can  be  tolerated  in  the  passive  mode 

This  could  be  explained  by  the  fact  that  there  are  two  errors  in 
determining  the  position  of  the  array  when  in  the  active  mode:  one  at  the 
time  of  transmission  and  the  other  at  reception  of  the  echo. 

The  effects  of  weighting  (Hamming)  in  the  presence  of  the  deterministic 
placement  errors  are  presented  in  Figs.  15  and  16  for  f  =  300  Hz  and  100  Hz 
respectively.  In  both  cases  the  weighting  function  smooths  out  all 
distortions  in  the  beam  pattern.  It  is  noticeable  that  with  f  -  100  Hz  the 
weighting  also  improves  bearing  resolution. 

The  effect  of  signal-to-noise  ratio  on  active  synthetic  aperture  has  also 
been  investigated.  The  results  are  very  similar  to  those  obtained  for  the 
passive  mode  (see  Fig.  1)  in  which  side  lobes  of  the  beam  pattern  are 
increased  at  low  signal-to-noise  ratios. 

A  beam  pattern  with  acceptable  main- lobe/sidelobe  ratio  was  obtained  in  the 
active  mode  for  a  signal-to-noise  ratio  not  lower  than  3  dB. 
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FIG.  13 

ACTIVE  SAS  -  EFFECT  OF 
DETERMINISTIC  PLACEMENT 
ERRORS 


FIG .  14 

ACTIVE  SAS  -  EFFECT  OF 
DETERMINISTIC  PLACEMENT 
ERRORS 


FIG.  15 

ACTIVE  SAS  -  EFFECT  OF 
WEIGHTING  IN  PRESENCE  OF 
DETERMINISTIC  PLACEMENT 
ERRORS  (Hamming  weighting) 


20 


SACLANTCtN  SR-91 


CONCLUSIONS 


The  results  obtained  from  computer  simulations  of  the  syntnetic  aperture 
techniques  as  applied  to  passive  sonar,  indicate  that  the  following  points 
need  consideration  in  designing  a  passive  synthetic-aperture  sonar  system: 

-  Signal  coherence  time  must  be  long,  or  at  least  equal  to  the 
total  synthesizing  time  necessary  to  form  a  synthetic  aperture. 

Signal-to-noise  ratio  at  the  output  of  the  synthesizing  array 
must  be  larger  than  3  dB  (see  para  1.5.2). 

'Deterministic'  placement  errors  caused  by  the  movement  of  the 
towing  vessel  allow  acceptable  beamforming  performance  if  the 
take-off  angle  of  the  track  is  not  larger  than  2°  and  the  total 
synthesizing  time  is  at  last  five  times  longer  than  the  period 
of  'oscillation'  of  the  vessel  around  its  nominal  straight 
course. 

-  Random  placement  errors  caused  by  the  fluctuations  in  the  motion 
of  the  synthesizing  array  can  be  tolerated  if  the  standard 
errors  are  less  than  0.16  . 

-  Weighting  (Hamming)  can  reduce  the  effect  of  placement  errors  on 
the  synthesized  beam  pattern. 

The  simulation  of  active  synthetic  aperture  indicates  that  the  following 
points  must  be  considered  in  designing  an  active  synthetic-aperture 
sonar  system: 

-  Maximum  operating  range  is  limited  but  can  be  increased  by 

lowering  the  frequency.  Moderate  ranges  of  10  to  15  km  can  be  1 

achieved  at  a  frequency  of  the  order  of  100  Hz.  Lower 

frequencies  may  require  long  synthesizing  arrays.  Higher 

frequencies  may  cause  diffraction  lobes  in  the  synthesized  beam 

pattern  at  the  ranges  considered. 

-  As  in  the  passive  mode,  a  signal-to-noise  ratio  of  3  dB  at  the 
output  of  the  physical  array  appears  adequate  for  beamforming. 

-  The  effects  of  placement  errors  is  more  pronounced  than  in  the 
passive  mode. 

-  Weighting  (Hamming)  minimizes  the  effects  of  placement  errors. 

-  The  angular  reflections  of  scatterers  must  be  coherent  during 
the  synthesizing  time. 
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These  results  indicate  that  the  coherence  time  of  the  signal  is  likely  to 
be  the  major  problem  in  the  passive  mode  if  large  synthetic  apertures  are 
required  at  low  speeds.  For  a  fixed  synthetic-aperture  length  an  increase 
in  speed  reduces  the  total  synthesizing  time  and  consequently  a  shorter 
coherence  time  is  required.  We  suggest  that  passive  sonar  should  make  more 
use  of  the  temporal  coherence  of  a  signal.  The  synthetic-aperture 
technique  looks  for  long  coherence  times;  perhaps  the  first  step  is  to 
seek  small  improvements  in  aperture  length. 

The  main  requirement  for  the  good  performance  of  an  active  synthetic- 
aperture  system  is  that  the  angular  reflection  characteristics  of  the 
scatterers  remain  coherent  throughout  the  synthesizing  time. 

In  the  passive  mode  there  is  no  limit  on  operating  rate,  because  the 
sampling  time  is  not  range  dependent  as  it  is  in  the  active  case.  The 
operating  range  for  an  active  synthetic  aperture  can  be  extended  by 
lowering  the  frequency. 

Low  signal-to-noise  ratios,  as  modelled  in  this  work,  have  little  effect  on 
the  bearing  resolution  achievable,  either  with  passive  or  active  synthetic- 
aperture  systems.  A  3  (IB  signal-to-noise  ratio  at  the  input  of  the 
beamformer  appears  adequate  for  both  passive  and  active  systems. 

The  random  placement  errors  can  be  made  negligible  by  operating  at  low 
frequencies.  The  'deterministic'  errors  can  considerably  affect  the 
performance  of  a  synthetic-aperture  system  if  the  take-off  angles  of  the 
moving  array  are  larger  than  2°  in  the  passive  mode  and  larger  than  0. b°  in 
the  active. 

One  very  important  conclusion  of  this  study  is  that  there  is  no 
experimental  information  in  the  open  literature  on  which  to  base  an 
accurate  evaluation  of  the  performance  that  could  be  achieved  with  an 
actual  system.  Measurements  are  required  that  are  directly  applicable  to 
the  synthetic-aperture  technique.  Consequently  they  must  be  made  with  a 
moving  receiver  platform  and  at  various  sampling  times.  Ultimately  the 
movement  of  the  source  must  be  included. 
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APPENDIX  A 

MATHEMATICAL  DESCRIPTION  OF  THE  PASSIVE 
SYNTHETIC-APERTURE  BEAMFORMING  PROCESS 


In  the  passive  synthetic-aperture  beamforming  process,  samples  of  the 
signal  received  at  successive  positions  of  the  synthesizing  array  as  it 
moves  through  the  synthesizing  sequence  are  compensated  for  the  delay  time 
between  samples  and  the  movement  of  the  array  and  then  added  coherently. 
The  prediction  of  performance  is  obviously  very  dependent  on  the  following 
assumptions: 

a.  The  radiating  source  is  stationary  and  radiates  with  a  constant 
amplitude  and  phase. 

b.  Own-doppler  nullification  (ODN)  is  carried  out  on  the  received 
signals  before  commencing  the  synthetic-aperture  beamforming 
process. 

We  start  by  analyzing  the  idealized  case,  in  which  it  is  assumed  that  the 
signa1  field  is  perfectly  coherent  in  space  ami  time  during  the 
synthesizing  period  (T  )  and  that  the  synthesizing  irray  moves  along  a 

perfectly  straight  course  at  a  known  constant  speed.  Ihereafter  we  examine 
the  effect  of: 

a.  lhe  temporal  signal  coherence  not  being  perfect  but  being 
characterized  by  a  coherence  time  (Tc). 

b.  Random  phase  and  amplitude  errors  in  the  signals  received  at 
successive  positions  of  the  synthesizing  array. 

c.  Random  and  deterministic  errors  in  the  position  of  the 

synthesizing  array  as  it  moves  through  the  synthesizing  sequence. 


a.i  idealized  BEAMFORMING 

On  the  basis  of  the  assumptions  given  above,  the  geometry  of  the  synthe¬ 
sizing  sequence  is  as  shown  in  Fig.  Al.  The  delay  and  steering  vectors 
required  to  focus  the  successive  signals  to  F  (the  assumed  position  of  the 
source)  can  be  derived  by  first  considering  the  received  phase  of  a  signal 
radiating  from  F  and  then  reversing  this  phase.  The  phase  of  the  signal 
received  at  the  start  of  the  synthesizing  sequence  is  taken  as  a  reference. 

Using  the  geometry  of  Fig.  Al,  the  phase  of  the  k^  signal  received  from  F 
is: 


4-Fk  =  u>|t1-(k-i)At  +  (KFk  '  Rfl)/c)  ,  (Eq.  A.l) 

where  wt^  is  the  reference  phase.  Consequently  the  delay  and  steering 
vectors  are  given  by 

dk  =  Ak  expl-juijt^Ck-Dit  +  (Rfk-RF1)/c)|  ,  (Eq.  A. 2) 

where  Ak  is  an  amplitude-weighting  function  that  can  be  used  to  shape  the 
synthesized  beam  pattern  if  required. 


Z5 


SACLANTCEN  SR-91 


Idealized  course  of  the 
synthesizing  array 

t 


FIG.  A1  GEOMETRY  FOR  PASSIVE  SYNTHETIC  APERTURE  MODE 

Initial  values  of  R  and  R  are  R  and  R  respectively . 
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MATHEMATICAL  DESCRIPTION  OF  THE  PASSIVE 
SYNTHETIC-APERTURE  BEAMFORMING  PROCESS 


In  the  passive  synthetic-aperture  beamforming  process,  samples  of  the 
signal  received  at  successive  positions  of  the  synthesizing  array  as  it 
moves  through  the  synthesizing  sequence  are  compensated  for  the  delay  time 
between  samples  and  the  movement  of  the  array  and  then  added  coherently. 
The  prediction  of  performance  is  obviously  very  dependent  on  the  following 
assumptions: 


a. 

The  radiating  source  is  stationary  and  ra 
amplitude  and  phase. 
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A.  1  IDEALIZED  BEAMFORMING 

On  the  basis  of  the  assumptions  given  above,  the  geometry  of  the  synthe¬ 
sizing  sequence  is  as  shown  in  Fig.  AI.  The  delay  and  steering  vectors 
required  to  focus  the  successive  signals  to  F  (the  assumed  position  of  the 
source)  can  be  derived  by  first  considering  the  received  phase  of  a  signal 
radiating  from  F  and  then  reversing  this  phase.  The  phase  of  the  signal 
received  at  the  start  of  the  synthesizing  sequence  is  taken  as  a  reference. 

Using  the  geometry  of  Fig.  Al,  the  phase  of  the  kth  signal  received  from  F 
is: 


<t,Fk  ~  +  (RFk  "  rfi)/c!  >  (Eq-  A1> 

where  iut^  is  the  reference  phase.  Consequently  the  delay  and  steering 
vectors  are  given  by 

dk  =  Akexp[-ju.{t1-(k-l)At*(RFk-RF1)/c}J  ,  (Eq.  A. 2) 

where  A^  is  an  amplitude-weighting  function  that  cm  lie  used  to  shape  the 
synthesized  beam  pattern  if  required. 
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The  next  step  is  to  define  the  signals  received  from  a  source  not  posi¬ 
tioned  at  F  to  which  the  above  vectors  will  be  applied.  By  varying  the 
bearing  of  this  source  it  will  be  possible  to  trace  out  the  synthesized 
beampattern. 

From  Eq.  A.  1,  the  signal  sample  obtained  from  a  source  at  position  S 
(see  Fig.  Al)  is  given  by: 

sk  =  Bk  exp[ jtu{tj-(k-l)At  +  (R$k~Rsi)/ci  1  •  (Eq.  A-3> 

In  this  expression  Bk  is  introduced  to  allow  for  the  signal  amplitude 

(bk)  and  for  the  amplitude  response  of  the  receiving  array  in  the 

direction  of  S,  and  to  normalize  the  amplitude  levels  relative  to  the 
amplitude  received  when  s  lies  on  the  broadside  axis  of  the  synthetic 
aperture.  Thus: 

Bk  =  bk  f(«s)  ■  f(Rs)  , 


where  f(a<.)  is  the  amplitude  response  pattern  of  the  physical  array  (a 

function  of  the  number  of  hydrophones,  spacing,  and  shading  function  of  the 
physical  array),  and  f(R<0  is  RSc/RSk'  Th’s  'atter  factor  is 

consequently  I  when  the  source  is  on  the  broadside  axis;  otherwise  it  is 
less  than  1. 


As  stated  previously,  to  derive  Eq.  A.  3,  we  assume  that  own-doppler 
nullification  has  been  carried  out  before  the  synthetic  beamforming  process 
starts.  However,  the  own-doppler  nullification  will  be  correct  only  for 
signals  from  the  focus,  since  this  is  the  direction  from  which  it  is 
assumed  the  signals  are  arriving.  Thus  the  present  analysis  makes  no 
allowance  for  the  actual  doppler  differing  from  the  assumed  value  as  the 
position  of  the  source  is  varied  to  generate  the  synthesized  beam  pattern. 

To  form  a  synthesized  beam  we  sum  N  samples  of  the  received  signals,  with 
the  appropriate  delay  and  steering  vectors  being  applied.  Using  Eqs.  A. 2 
and  A. 3  the  beam  output  voltage  is  given  by: 


N  N 

V(P)  =  2  Sk  •  d  k=  2  S 


k=l 


k=l 


k  ' 


Thus 


V(P)  =1  Ak  Bk  exp(  j(0  Sk-0fk) ) 


where 


(Eq.  A. 4) 


(Eq.  A. 5) 


Sk  "  (RSk'RSlWc 

(Eq. 

A. 6) 

Fk  T  <RFk'RFl  )uj/c  • 

(Eq. 

A.7) 
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Then  in  the  usual  manner  we  can  write  for  the  beam  output  power: 

N  N 

P  t } V -V*}  -  t { |  i  Ak  Bkcos(0sk-UF|<)|2  M  i.  A  kB  ksin(uSk-MFk)|2| 

k=  1  k-  1 

Since  this  contains  no  time-varying  components  we  have: 

P=  AkBkcos(e5k-0Fk)l2Mi  A.S.sin^-ep,  )l2  • 

k=1  k=1  (Eg.  A. 8) 


Referring  back  to  Eq.  5  it  can  be  seen  from  Fig.  A. 1  that: 

Rpl  =  [ (L/2  +  Rq  tan  pp)2  +  RQ2]h  (Eq.  A. 9) 

Rpk  =  [{L/2  +  Rq  tan  pp  -  (k-l)UAt!2  +  Rq2]^  (Eq.  A. 10) 


and  that  the  equations  for  Rc;l  and  R,..  are  exactly  the  same,  with  pF  being 
replaced  by  (S<.. 

Returning  for  a  moment  to  the  equation  for  V(p),  we  have  from  Eqs.  A.E>,  A. 6 
and  A. 7  that 

N 

V(fi)  =  ^  A  kB  kexp[j{RSk-  R  S1+ Rpi  -  Rp^iu/cl  (Eq.  A. 11) 

Then  using  Eqs.  A. 9  and  A.  10,  and  the  similar  ones  for  R<-  and  Rsk  it  can 

be  shown  after  some  tedious  algebra  that  for  R  >>L,  namely  in  the  far-field 
of  the  synthesized  aperture:  0 


V(p)  =  exp[ j  ^  (sin  fSp  -  sin  fc)]  •  H  ^  . 

k=l 

•  exp[j^  (k-1)  UAt(sin  pp  -  sin  p$)] 
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This  is  the  standard  expression  for  the  far-field  beam  pattern  of  a  linear 
array  of  N  hydrophones  equally  spaced  at  intervals  of  Uat. 

It  is  usual  to  normalize  P  by  dividing  it  by  N2 ,  so  that  its  maximum  value 
is  1  when  {3p  equals  pk-  However,  for  this  present  study  the  normali¬ 
zation  has  been  done  with  respect  to  N  only.  Thus  with  =  1  for 

all  values  of  k,  P  will  have  a  maximum  value  of  10  log  N  when  Pp  equals 

h- 


A. 2  EFFECT  OF  THE  TEMPORAL  COHERENCE  OF  THE  RECEIVED  SIGNAL 

Equations  A.  3  and  A.  8  for  the  kth  signal  sample  and  the  synthesized  beam 
power  respectively  are  both  based  on  the  assumption  that  all  the  received 
signal  samples,  after  phase  correction  for  movement  and  time  difference, 
are  fully  coherent.  This  means  that  the  normalized  temporal  correlation, 
or  coherence  in  our  case  of  a  narrowband  signal,  is  1  during  the  period  T  . 
Namely,  that: 

E{Sk  ’  Sk+i! 

- - - =  p(S.  )  =  1  ,  (Eq.  A. 13) 

[E{S2}  •  E{S2k+ii ]”* 

where  S^,  Sk+^  are  defined  by  Eqs.  A. 4,  A. 5,  A. 6  and  A. 7. 

In  order  to  examine  the  likely  effect  of  a  lack  of  coherence,  we  assume 
that  there  is  an  exponential  decay  with  time  of  p(st).  as  follows: 

p(St)  =p(St)  •  exp[-|tk+1  '  t  j/T  ]  ,  (Eq.  A. 14) 

actual 


where  Tc  is  referred  to  as  the  coherence  time  of  the  signal. 

In  other  words  we  have: 

E{Sk'Sk+i*  =  VSk+i  exp^"  ltk+i‘tkl/V  ’  rEq-  A'15) 


The  equation  for  the  synthesized  beam  output  power  is  now: 

N  N 

E{(  2  S  )  •  (  2  S.  )*)  , 

k=l  K  k=l  K 


thus: 


P  = 


N  N-l  N-i 

2  S2  +  22  2 

k=l  K  i=l  k=l 


E{S. 


n+i* 


N  N-l 


=  2  S2  +  2  2 
k=l  K  i=l 


N-i 

2 

k=l 


Sk  Sk+i 


exp[-|tk+i 


l/TcJ- 


(Eq.  A. 16) 
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A.  3  EFFECT  OF  RANDOM  PHASE  AND  AMPLITUDE  ERRORS  IN  THE  RECEIVED  SIGNAL 

Random  fluctuations  can  occur  in  the  amplitude  and  phase  of  the  received 
signal  due  to  a  number  of  fac1  rs.  Examples  of  some  of  these  factors  are 
fluctuations  of  the  radiating  source,  fluctuations  caused  during  the 
propagation  of  the  signal,  the  presence  of  noise,  and  multipaths  causing 
constructive  and  destructive  interference  along  the  synthesizing  path. 

It  is  not  possible  to  characterize  these  fluctuations  from  the  published 
information  on  fluctuations,  because  none  of  the  information  is  strictly 
applicable  to  the  synthetic-aperture  technique.  All  the  published 
information  has  been  obtained  with  stationary  receivers  and  some 
measurements  were  made  with  spaced  receivers.  Flowever,  these  were 
simultaneous  measurements  and  therefore  not  useful  for  synthetic-aperture 
studies. 

Consequently,  we  base  the  following  analysis  on  the  assumption  that  such 
random  errors  can  be  studied  by  including  random  fluctuation  components,  q 
and  c  ,  in  the  amplitude  and  phases  respectively  of  the  received  signal. 
We  further  assume  that  q  and  e  are  independent  of  each  other  and  of  the 
actual  amplitudes  and  phases,  and  can  be  characterized  as  first-order 
Markov  processes.  Then  in  order  to  quantify  the  effect  of  these 
fluctuations  (bearing  in  mind  the  lack  of  experimental  data)  their  root 
mean  square  (rms)  values  have  been  related  to  the  signal/noise  ratio.  The 
complete  assumptions  can  be  written  as: 


q£  -  q2k+i  noise  power 

Eink  •  nk+i!  =7^kexP[-|tk+.  -  \ J  /te  1 

IT  ^k+i 

T*  =  0  __ 

E{tk  *  ek+i}  =  ‘V^I'K+i  •  tkl/TEl 


(Eq.  A. 17) 
(Eq.  A. 18) 

(Eq.  A. 19) 

(Eq.  A. 20) 


where 

S/N  =  signal-to-noise  power  ratio  at  the  output  of  the  synthe¬ 
sizing  array 

=  correlation  time  for  the  amplitude  and  phase  fluctuations, 
which  are  assumed  to  be  equal. 


As  stated  earlier,  the  need  for  the  above  assumptions  highlights  an  area  of 
considerable  weakness  in  our  knowledge  of  acoustic  reception  in  the  sea 
from  a  moving  platform. 
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Using  the  above  assumptions  we  can  first  rewrite  Eq.  A. 8  fot  the  power 
output  of  the  synthesized  beam  as: 

N 

P  =  E{[^2  Ak  f(ok  )  f(R  k  )  •  (bk+nk>  cos(8k+ek)]2l  + 

N 

+  E{[  ak  f(ak)  f(Rk  )  •  (bk+nk)  sin(6k+ek)]2}.  (Eq.  A. 21) 

k=l 

This  can  be  put  in  the  form  of: 

N  N 

P  =  E{[  2.  r.  ]2}  +  E{[  2  pk)2i 
k=l  K  k=l 

where 

rk  =  Ak  ‘  (W  *  * 

pk  =  Ak  f(“k)  f(Rk5  ‘  (bk+nk)  ‘  sin(6k+£k) 

and 

8k  =  0^k  -  6j-k  .  (see  Eqs.  A.  5  to  A.  8) 

Equation  A. 21  can  now  be  restructured  to  read: 

N  N-l  N-i 

P=2  [E{r2}  ♦  E{p2J]  +  2  i  2  [E(rk-rk+i}  h.  E{pk-pk+.}]  .  (Eq.  A. 25) 

k=l  i=l  k=l 

From  Eqs.  A. 23  and  A. 24  the  first  summation  in  Eq.  A. 25  becomes 

N 

^EUAkf(ok)f(Rk)  (bk*nk)l2!  . 

which  we  will  write  as: 

N 

1  q2k  ' 
k=l  * 

Similarly, 

E{VW  =  VVi  lEfcos«YCk)  COs(eka+tk+i)!1  (Eq-  A.  28) 


(Eq.  A. 26) 


(Eq.  A. 27) 


(Eq.  A. 22) 

(Eq.  A. 23) 
(Eq.  A. 24) 
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Thus  finally  we  have: 

N  _ 

P  =  i  A*  f2(ak)f2(Rk)(b2k+n2k)  ♦ 

+  2i=\  ^'Akf(“k)f(Rk)-Ak+if(“k+i)  f(Rk+i>  ’lbk  + 

+  ^k  exp('l  1  k+i-tkl/T E> 1  ‘ 1  6XP('I  1  k+i_t  kl/T  E  )  1  COS(0  k+0  k+i)} 

(Eq.  A. 34) 


A.  4  EFFECT  OF  PLACEMENT  ERRORS 

The  expression  "placement  error"  is  used  to  describe  the  difference  in 
position  between  the  actual  and  assumed  positions  of  the  synthesizing  array 
as  it  moves  through  the  synthesizing  sequence.  Such  errors  will  cause  the 
delay  vectors  to  be  incorrect.  Thus,  for  example,  a  source  at  the  "focus" 
of  a  synthesized  aperture  will  not  be  focussed  correctly.  It  should  be 
noted,  however,  that  since  the  degradation  caused  by  these  errors  is 
obviously  dependent  on  the  ratio  of  the  placement  error  to  the  wavelength 
in  use,  their  effect  can  be  reduced  by  reducing  the  operating  frequency. 

The  present  analysis  attempts  to  simulate  two  types  of  placement  error; 
Firstly,  a  purely  random  positional  error,  e.g.  random  navigational  errors 
or  an  erratic  movement  of  the  synthesizing  array  through  the  sea  and, 
secondly,  a  placement  error  caused  by  a  non-linear  track  resulting  from  the 
oscillations  of  a  moving  platform  about  the  desired  straight  course  at  sea. 
The  second  type  of  error  is  considered  deterministic  because  the  actual 
value  of  each  successive  error  can  be  mathematically  defined  whereas  the 
first  type  must  be  treated  statistically. 


A. 4.1  The  Effect  of  Random  Placement  Errors 

The  purely  random  placement  errors  are  modelled  by  assuming  a  gaussian 
distribution  of  actual  array  positions  about  the  assumed  array  position, 
with  the  errors  being  independent  of  each  other  and  independent  of  the 
actual  positions.  The  gaussian  distribution  is  assumed  to  be  circular,  so 
that  it  produces  the  same  effect  for  all  signal  directions.  Thus  it  is 
defined  by  a  single  variance  (o2)  that  can  be  converted  to  a  variance  in 

phase  by: 

?n  2 

VAR  (Phase)  =  (o  •  ,  )  (Eq.  A. 35) 

y  a 
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Then,  using  £q.  A.  34  and  the  assumptions  outlined  above,  we  obtain  the 
analytical  expression  for  the  synthesized  beam  power  output  as: 


P  = 


N 

Z  A2  B2 
k=  l  k 


N-l  N-i 
2  Z  Z  [A 
i=l  k=l 


Bk  Ak+i  Bk+i  '  (1_VAR(Phase))  c°s(0 k~ek+i ) ]  . 

(Eq.  A. 36) 


A . 4 . 2  Effect  of  Deterministic  Placement  Errors 

As  stated  earlier,  the  aim  is  to  simulate  the  type  of  non-linear 
oscillating  track  that  usually  occurs  at  sea  when  a  moving  platform 
attempts  to  follow  a  straight-line  course. 

We  therefore  assume  that  the  speed  of  the  synthesizing  array  along  the 
oscillating  track  is  constant  but  not  necessarily  the  assumed  speed,  and 
that  the  oscillation  can  be  characterized  by  an  across-the-track  speed 
component  of: 

U  (t)  -  U  sin  y  ■  cos(2nt/TQ)  ,  (Eq.  A. 37) 

where 

tl  =  actual  speed  (constant)  along  the  track 

y  =  maximum  off-track  angle  of  the  oscillating  track. 

Tq  =  period  of  the  oscillation  (see  Fig.  Al) 

Thus  for  the  along-track  component  of  velocity  we  have: 

Ux(t)  =  |U2  -  U 2 ( t ) ] 15  =  till  -  cos2(2nt/To)]1'5  (Eq.  A.  38) 

using  the  assumption  that  sin2y  •  cos2(2nt/TQ)  is  always  <<1. 

Consequently  the  deterministic  placement  errors  are  given  by  the 
displacement  of  the  oscillating  track  from  the  nominal  straight  line,  as 
fol lows: 


Tk 

AX.  =  X  0  (t)  •  dt 
TOS  X 


(Eq.  A. 39) 
Tk 

Ay.  =  X  U  (t)  •  dt 
K  TOS  y 


where 

Tk  = 

time  corresponding 
synthesizing  array 

to 

the  successive 

positions  of  the 

TOS  = 

time  corresponding 

to 

the  initial 

position  of  the 

synthesizing  array  on 

the 

oscil lating  course. 

34 


SACLANTCEN  SR'91 


The  effect  of  these  deterministic  errors  is  then  introduced  by  modifying  in 
Eq.  A.  11  the  values  for  the  ranges  R,,  and  R$.  (see  Eqs.  A. 9  and  A.  10). 
Thus  R^  and  R^k  become: 

R'sl  =  [(L/2  +  Rq  tan  ps)2  +  (Rq  +  Ayp2]'"5  (Eq.  A. 40) 

R'k=  ( (L/2  +  Rq  tan  ps~  Axk)2  +  (RQ  +  Ayk  )2  ]H  (Eq.  A.  41) 


Therefore  in  the  Eq.  A. 34  of  P  we  will  have 


6,1  =  6'  ,  -  8, , 

k  Sk  Fk 


(Eq.  A. 42) 


6k+i  =  ®  Sk+i  "  6  Fk 


(Eq.  A. 43) 


where  6'sk  and  9'jk+^  are  obtained  from  Eqs.  A. 6  and  A. 7  using  R'S1 

R'  .  So  that  with  deterministic  errors  alone  Eq.  A. 34  can  be  written: 
Sk 


and 
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Idealized  course  of  the 
synthesizing  array 


El  Transmission 

Expected  reception  period  of  echo  from  F 
RiSil  Actual  reception  period  from  scatterer  S 

FIG.  D2  TIMING  SEQUENCE  OF  ACTIVE  IDEAL  BEAMFORMING 
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The  k-th  echo  received  from  the  focus  F  is  given  by: 


sk  =ak'  exp|ju)[(t1+  1  (At).  F  +  (k-l)t)  +  (R'kF+  R"kF  )/c] 
I  k=l 


s.  =  a  •  expi ju(t  +T  p)i 


(Eq.  B. 4) 

(Eq.  B. 5) 


To  bring  these  signals  into  phase  with  the  first  received  signal  thus 
requires  a  steering  vector  of  the  form: 


Dk  =  exp(-jujTkF)  , 


where 


T 


kF 


N 

I  (At)  .  p  +  (k-l)t 
k=l 


(Eq.  B. 6) 

(Eq.  B. 7) 


The  actual  phase  shift  of  the  echo  from  the  scatterer  S  relative  to  the 
first  transmission  (t  =  t^)  is  given  by: 

rk  =  b  k'  exPiJl°(ti  +  TkS  •  (Eq-  B.8) 

where  (see  Fig.  B2) 

TkF  =  TkS  *  (AtkF  ‘AtkS)  ‘  (Eq-  B  9) 


To  form  beams  we  require: 

N  N 

v=  2_^rk  ‘  Dk=  ^  b  kexP{  ju»Tks  }  ak  exp{-jwTkF}  (Eq.  B.10) 

N 

=  i  \  •  ak  •  exp{juo(TKS-Tkp)}  (Eq.  B.ll) 

N 

=  2  w  k  exp{  je  k]  ,  (Eq.  B.  12) 

where 

bk  =  amplitude  response  of  the  receiving  physical  array  in  the 

direction  of  the  scatterer  S. 

ak  =  amplitude  shading  factor  that  can  be  used  to  suppress  the 

side-lobes  of  the  synthesized  beam,  if  required. 
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T 


“TkF  = 

phase  shift 
to  focus  the 

of  the  echo  relative  to  first  transmission,  required 
synthesized  beam  to  position  F. 

wTkS  = 

actual  phase 

shift  of  the  echo  from  the  scatterer  position  S. 

mk  = 

bk  *  ak 

(Eq.  B. 13) 

IU  = 

2nf 

(Eq.  B. 14) 

where  f  is  the  frequency  of  the  transmitted  signal. 


The  phase  shift  6^  of  Eq.  B.12  is  calculated  based  on  the  geometry  of 
Fig.  B1  and  the  timing  sequence  of  Fig.  B2. 

From  the  relations  (Fig.  81) 


(,AtkF)2  ^  R2kS-  2VAtkF  R‘‘s  •  cos<*/2  -  orkp)  =  R'* 


At 


RkF  +  RkF 


kF 


(Eq.  B. 15) 

(Eq.  B. 16) 


we  obtain 


R" 

kF 

V 

kF 


hi  -  -  *  -hfh;h 


‘•a) 


(Eq.  B. 17) 


From  Eqs.  B.15,  B.16  and  B.17  we  derive: 


“kF ■  {‘♦HI  -  •*»*)•  {[Hi •  •<»„)]*  * 

[>-<D2]  •  [hi  -  ^ ”»*)]♦  i}'-  Mi)2]'1} 

(Eq.  B. 18) 

From  Fig.  B1  we  also  derive: 

RkF  +2vTkF  RkF  si"“kF  +  <vV2=RZoF 

(Eq.  B. 19) 

RkF  =  *  vTkF  sinakF  +  {R2oF'  (vTkF  C0S“kF)2!!S 

(Eq.  B. 20) 
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The  angle  akp  is  calculated  from  (see  Fig.  Bl): 

R  tan  6  =  R  tan  8C  +  R  tan  a.  r  +  vT. c 
o  'm  o  ^F  o  kF  kF 

R  tan  B  =  v  •  T/2  , 

o  rm 

where  T  is  the  total  integration  time. 

We  obtain  for  u.  r: 

kF 


2T 


kF, 


okp  =  arctan  J(1 - — )  tan  Bm  -  tan  Ppj 


where  Tkp  is  calculated  by  Eq.  B. 7. 
For  Atk<.  we  obtain  similarly: 

At 


kS 


wi  th 


•[>*?(?-  -*'”•«)] 

[*"(')*]■['('  '  2  SI"“b)]  *  1  f  '  [‘'(O’]"1] 


RkS  =  '  vTkS  s1nukS+  lRV  (vTkF  cos“kS)2}!5 
and  the  angle  «kS  is  given  by: 

2TkF 

«kS  =arctani(l - p~  )  tan  Pm  -  tan  Ps) 

The  phase  shift  0,  is  calculated  from  Eq.  B14  using  Eqs.  B.9 
B. 24.  K 


Equation  B.20  can  also  be  written  for  the  generic  time  t  =  tk  in 


Rk  =  ‘  vtk  sin 


"  ,  /  (vtkcos  OzVl  ** 

°k  t°  \l — ^r~). 


In  the  approximation  of  far-field  (i.e.  Vtk/RQ  «  1)  we  have: 


Rk  ~  ”  vtk  s1n  uk  +  R 


(vtkcos  a k)2" 
o  1  '  W 

L  o 


(Eq.  B. 21) 
(Eq.  B. 22) 

(Eq.  B. 23) 

+ 

(Eq.  B. 24) 

(Eq.  B. 25) 

(Eq.  B. 26) 
,  B.18,  and 

the  form: 
(Eq.  B. 27) 

(Eq.  B. 28) 
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The  synthesized  beam  voltage  is  given  by  the  sum: 

v= irkexp[-^(Av  -  AtkS>] 

If  we  assume,  in  the  far-field,  ak<. 

(sin  akS  +  sin  akf:)  •  (sin  akS  -  sin  ak(;)  2  sin  akf:(sin  akS  -  sin  akf) 


(Eq.  B. 38) 


Therefore,  from  Eqs.  B.37  and  B.38  we  obtain: 
V  ! 


wk  ex?  '  ^  TkF  (l  +  T  Sin  “  kf)  •  (  Si  1 "  “  kS  '  sin  “  Kf) 


Under  the  assumption  of  ^  sin  akp  <<  1 


V  “^expJ-^T  kF-  (sinakS-  sin^)] 


(Eq.  B. 39) 


(Eq.  B. 40) 


(Eq.  B. 41) 


where  Tkfr  is  given  by  Eq.  B.7. 

In  the  far-field,  the  difference  between  the  sampling  intervals  becomes 
very  small.  Therefore,  if  we  assume  that: 


TkF  =  kAt 

and  in  the  small  angle  approximation: 


* 


kS  2 


T  -  P 


m  N  HS 


’kF  2  *  P  m‘  V  PF 


(Eq.  B. 42) 

(Eq.  B. 43) 

(Eq.  B. 44) 


we  obtain  for  the  beam  voltage: 

N  n  2  i 

V  -  ^exp  kAt  *  (Sln  “kS  ’  s1n  akF ^ 


(Eq.  B. 45) 
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The  synthesized  beam  pattern  is  then: 


V-V*  = 


.  /Nun/Ar  \ 

s  l  n  I  — —  s  l  n  t|i  J 

s  in  ^  sin  t(i^ 

.  ujvAt  •  ,  "X 

= 

.  (nd  .  ,3 

sin  “E —  sin 

sin^—  sin  till 

(Eq.  B. 46) 


which  is  the  standard  form  for  a  linear  array  in  far-field  conditions  and 
composed  by  N  hydrophones  uniformly  spaced  by  d. 


Equation  B.46  represents  the  beam  pattern  formed  with  2N  samples.  That  is, 
the  effective  synthesized  aperture  is  2Nd.  Consequently  the  active 
synthesized  beam  pattern  has  a  3  dB  beamwidth,  which  is  half  of  that  for  a 
normal  passive  array  of  the  same  length.  Also,  as  discussed  in  par.  2.4, 
the  active  synthetic  aperture  array  differs  from  the  normal  active  array, 
which  would  have  a  beamwidth  forty  percent  greater  for  the  same  length. 


B2  Beamforming  in  the  presence  of  random  amplitude  and  phase  errors 


The  synthetic  beam  former  makes  the  sum  of  N  echoes: 

N 

V  =  2  s.  (Eq.  B. 47) 

k=l 


The  k-th  component  is  given  by: 

rk  =  wk  exp(j0k)  ,  (Eq.  B.48) 

with: 

wk  =  WR  +  nk  .  (Eq.  B. 49) 

where  nk  is  a  random  amplitude  fluctuation  term. 

And  with: 

6k  =  6k  +  ek  (Eq.  B‘ 

where  ek  is  a  random  phase  jitter. 

From  Eqs.  B.47,  B.48,  B.49  and  8.50  we  have: 

N 

V  -l  (wR  *  nk)  exp(  j(0  k  +  e  k  ))  (Eq.  B.51) 
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rk  =  <  W  cos(V£k} 


(Eq.  B. 58) 


SACLANTCEN  SR-91 


and 

Pk =  (W  S1n(ek+Ek) 


(Eq.  B. 59) 


For  the  analysis  of  beamforming  we  assume  that  the  two  random  variables 
and  rjk  are: 

a.  Mutually  independent 

b.  Correlated  in  time 


Under  the  first  assumption,  from  Eqs.  B.58  and  B.59  we  have: 


S_ 

UJ 

=  E{  (w 

k+,l 

k>2* 

•  E{cos2(0k+ck)j 

(Eq. 

B.  60) 

E^ 

=  E{  (w 

k+P 

k>2> 

•  E{sin2(ek+tk)j 

(Eq. 

B.  61) 

and 

E*Wi} 

=  E{(w 

k+P 

k)  cos(W  ‘  (Vi+\H) 

cos(6 

k+i+Vi)} 

(Eq. 

B.  62) 

E{pk-Pk>i} 

=  E{  (w 

k+P 

k>  s: 

in(W  ‘  (Vi+Vi} 

sin(6 

k+i+ek+i^ 

(Eq. 

B.  63) 

or 

Efrk'Vij 

=  E{(w 

k+P 

k>  ' 

(WVi)!  '  E{C0S(6 

k+£k> 

•  cos(W 

Vi  » 

(Eq. 

B.  64) 

E^k'Pk+i} 

=  E{(w 

k+rl 

k>  ■ 

(WVi)]  *  E{sin(e 

k+£k> 

•  sin<W 

Vi  )! 

(Eq.  B. 65) 


The  expected  beam  power  output  is  therefore  given  by: 

E{ V* V*}  :  iqy  zV  I=jqkqk+i  •  [Eicos^e,)  •  cosO^e^)}  ♦ 


E{sin(ek+ek)  •  sin(^+j+Vj)}] 

where 

<»k  «  E{(wk+n^}  ;  qkqk+i=  Ef(wk+nk)'(wk+.+nk+j)} 


(Eq.  B. 66) 

(Eq.  B. 67) 
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Considering  the  trigonometric  relation: 

cos(0k+ek)  •  cos(6k+i+fW  +  sin(ek+£k)  •  sin(ek+i+ek+i} 

-'cos(ek+£k‘ek+r£k+i)  C0S|(V^i  )+(£k'Vi )] 

,c°s(ek-0k+j)  •  cos(£jc  -£k+i  )  -  sin<Vek+i>  sin(ek’tk+i) 

(Eq.  B. 68) 

Equation  B.66  can  be  written: 

EiV-V*}  =  i  d*  2V  2  qkQktj-  l««<Vek+  ,  )  •  E{cos(£  k-e  k+i)i  - 
k=l  i=l  k=l 


sin(Ok-0k+i)  •  E(sin(VVi  )}] 


(Eq 

.  B. 69) 

In  the  approximation  of 

small  errors: 

E{sin(ek-tk+1)i  ^ 

E!ek 

ek+ii=  0 

(Eq 

.  B. 70) 

E{cos(£k'ekri)^ 

EU  - 

(Vck*i>\ 

2  '  ' 

!*? 

IE‘eV  +  E<eW 

-  2Eiek-e 

k+i>) 

(Eq 

.  B. 71) 

The  expected  value  of  the  beam  power  output  is: 

N  N-l  N-i 

E!V.V*i  =1^2.^  ^qkqk+1.{cos(ek-8k+i)  - 

[1-  |  E{r.2j  -  |  E{e2k+i}  +  E{ek-ek+i}])  (Eq.  B.72) 
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It  is  assumed  that  the  random  phase  jitter  ck  and  the  random  fluctuation  in 
the  amplitude  qk  are  described  by  a  stationary  Markov  process  of  the 
first  order.  Thus  the  conditions  on  and  ok  are: 

Zero-mean 


Eiek!  =  E!nk!  =  0 

(Eq.  B. 74) 

with  variances: 

=  a2£k  ■  E^  s  °2nk 

(Eq.  B. 75) 

and  correlation  given  by: 

E*VEk+i*  =  ae.  ’  e*P[-IVW/Tll 

k 

(Eq.  B. 76) 

E{Vnk+i!  =  °nk  ‘  exp[-|tk-tk+.|/t2l 

(Eq.  B. 77) 

with  t ^  and  time  constant  of  the  fluctuations  of  phase  and  amplitude 

respectively. 

Stationarity  of  phase  and  amplitude  fluctuations  implies  that: 


E{^kl  = 

=  E!^k+ij 

(Eq. 

B.  78) 

and 

EK»  = 

:  E^W 

(Eq. 

B.  79) 

In  the  condition 

of  a  stat 

ionary 

signal ,  i . e. 

EK>  = 

:  EK+1-» 

» 

(Eq. 

B.  80) 

by  combi n 

ing  Eqs 

.  B.  74  to 

B.  80 

we  obtain  for 

the  expected  beam 

power 

output 

the  expression: 

N 

N- 

1  N-i 

E{ V-V* j  = 

*(E{w 

k=l 

-2j  +  a2  ) 

k  V 

+  22 
i= 

2  [Eiw2i  + 
1  k=l 

°nkexpH  VVi| 

/\}] 

•u  - 

a2  (1  - 

exp{-| 

tk‘tk+i|  /t2})1 

K 

1 

(Eq. 

B.  81) 

where  0k  is  the  phase  calculated  from  the  geometry  of  Fig.  B1  in  the 
absence  of  any  errors  (see  Eq.  B.12). 
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83  Analysis  of  Placement  Errors 

Based  on  the  geometry  of  Fig.  63  we  calculate  the  placement  error  given  by 
the  di fference: 

AkT  s  <r'kT+r"kT>-  <RkT  *RkT  >  B82> 


The  second  quantity  of  the  right  side  of  Eq.  B.82  is  defined  by  the 
geometry  of  the  ideal  case  of  straight-line  course  (see  Fig.  Bl). 

The  first  quantity  is  computed  from  (see  Fig.  B3): 


rkT  = 

s'i? 

+  R'kT'  26  k 

R‘kTC0SKs  +  \ 

+  arctan 

yk  1 

TTk  1 
xk 

(Eq.  B.83) 

rkT  = 

6"k 

+  R"kf  26"k 

rV0SKs  +  2 

+  arctan 

yk* 

x5*  ] 
xk 

(Eq.  B. 84) 

where 

“‘kS 

and 

a"kS  are 

the  aspect  angles  of  the  target  as 

defined  in  the 

straight-line  path  geometry. 


6k  =  [y*k2  +  (Xk  '  Xk)2)!*  (Eq.  B.85) 

6n  s  [y**2  +  .  x**)2]!}  (Eq.  B.  86} 


y 


FIG.  B3  GEOMETRY  OF  PLACEMENT  ERRORS  (Active  operation) 
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